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ABSTRACT 



Context. Recently, unanticipated magnetic activity in ultracool dwarfs (UCDs, spectral classes later than M7) has emerged from a 
number of radio observations. The highly (up to 100%) circularly polarized nature and high brightness temperature of the emission 
have been interpreted as requiring an effective amplification mechanism of the high-frequency electromagnetic waves - the electron 
cyclotron maser instability (ECMI). 

Aims. We aim to understand the magnetic topology and the properties of the radio emitting region and associated plasmas in these 

ultracool dwarfs, interpreting the origin of radio pulses and their radiation mechanism. 

Methods. An active region model was built, based on the rotation of the UCD and the ECMI mechanism. 

Results. The high degree of variability in the brightness and the diverse profile of pulses can be interpreted in terms of a large-scale 
hot active region with extended magnetic structure existing in the magnetosphere of TVLM 513-46546. We suggest the time profile 
of the radio light curve is in the form of power law in the model. Combining the analysis of the data and our simulation, we can 
determine the loss-cone electrons have a density in the range of 1.25 x 10 5 - 5 x 10 5 cirT 3 and temperature between 10 7 and 5 x 10 7 
K. The active region has a size < 1 Rj up , while the pulses produced by the ECMI mechanism are from a much more compact region 
(e.g. ~0.007 ^j up ). A surface magnetic field strength of «7000 G is predicted. 

Conclusions. The active region model is applied to the radio emission from TVLM 513-46546, in which the ECMI mechanism 
is responsible for the radio bursts from the magnetic tubes and the rotation of the dwarf can modulate the integral of flux with 
respect to time. The radio emitting region consists of complicated substructures. With this model, we can determine the nature (e.g. 
size, temperature, density) of the radio emitting region and plasma. The magnetic topology can also be constrained. We compare 
our predicted X-ray flux with Chandra X-ray observation of TVLM 513-46546. Although the X-ray detection is only marginally 
significant, our predicted flux is significantly lower than the observed flux. Further multi-wavelength observations will help us better 
understand the magnetic field structure and plasma behavior on the ultracool dwarf. 

Key words, magnetic fields - radio continuum: stars - stars: low-mass, brown dwarfs - polarization - masers 



1. Introduction 

A large number of recent radio observations indicate that intense 
k^I ' magnetic activity exists in ultraco ol dwarfs (UCDs), i. e. objects 
. , with spectral type later than M7. iBerger et al.l d200lh reported 
■ the first detection of quiescent and flaring radio emission from 
?-h \ the M9 brown dw arf LP944-20, in wh ich a bright X-ray flare 
, 5^ i was also detected (Rutledge et al.ll2000h . with anomalous quies- 
cent radio luminosity at least four orders of magnitude larger 
than predicted from an empirical relation between the X-ray 
and radio luminosities of activ e stars with spectral types from 
F to M dGuedel & Benzlll993l) . The detection of electron cy- 
clotron maser (ECM) emission provided the first confirmatio n 
of kilogauss fields for a late M dwarf (iHallinan et al. 2006), 
and subsequently led to the discvery that even cooler L type 
dwarfs can also pos sess magnetic fields in the kilogauss range 
dHallinan et alJ2008h . Radio observations exclusively enable the 
measurement of magnetic fields on coo l, brown dwa rfs, and pos- 
sibly also for the very faint exoplanets (Zarka 2007). 

Chromospheric Ha emission and coronal X-ray emission 
show a sharp decline in Lna/Lboi and Lx/^boi beyond spectral 



type M7 dNeuhauser et all Il999b iGizis et al.l |2000t IWest et al.l 
l2004HSterzer et al.ll2006aHSchmidt et alJl2007h" which would be 
consistent with lo wer fractional ioniz ation in atmospheres of 
later spectral type (Mo hantv et al.ll2002l) . Recently, however, a 
lot of evidence, such as quiescent and flaring Ha emission from 
some L and T dwarfs (Reidetal. 1999; Burgass er et al. 2000; 
Lieb ert et al.l 12003b iRei ners & Basri 2007; Rockenfelle r et al .1 
boOrjllStelzer et all2006bH Schmidt et alJl2007l) . FeH fines from 
cool M dwarfs (lAfram et al.ll2 009). strong X-ray emission from 
one such source dAudard et al.ll2.007b , all suggest t hat perhaps an 
efficient magnetic dynamo could be operational (| P aikerlll955h 
in a fraction of the UCD population (IReiners & B asri 2010J) and 
that magnetic reconnection events could occur on these kinds 
of cool objects. Such magnetic activity and thus the strong ra- 
dio emission could be associated with the differential rotation 
between the atmosphere and the core of the UCD. 



iBergerl (l2002h reported Very Large Array (VLA) observa- 
tions of 12 late M and L dwarfs in the solar neighborhood. 
Flare-like outbursts as well as persistent quiescent emission 
were detected from three of the 12 sources, TVLM 513-46546 
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(TVLM 513 hereafter), 2MASS J00361617+1821 104 (2MASS 
J0036+18) and BRI 0021-0214. Among these three radio-active 
sources, TVLM 513 and BRI 0021, plus the first detected 
radio-active UCD, LP944-20, all have ra pid rotational velocities, 
v sin i >30 km s _1 . lHallinan et al.l (120061) presented observations 
of the rapidly rotating M9 dwarf TVLM 513 obtained simulta- 
neously at 4.88 and 8.44 GHz using the VLA. The periodic radio 
emission at both frequencies indicated a period of ~2 hr in ex- 
cellent agreement with the estimated period of rotation of the 
dwarf based on its v sin ; of ~60 k m s _1 . One more radio-active 
source from lHallinan et al.l d2008l) . LSR J1835+3259, also has a 
high rotational velocity, v sin i «50 km s _I . 

Up to now, about 10 radio active 
with spectral type from M8 to L3.5, 
nary system (2MASS J07464256+2000321, 
0746+20), have been found from various 



ultracool dwarfs, 
including a bi- 
hereafter 2MASS 
su rveys dBergerj 

2002; [Burgasser & Putmanl l2005t i Berged 120061: lHallinan et all 
2006, 120071 12008b lAntonova et alJl2008h . Three of these have 



been shown to have periodi c radio emissi on, with periods of 
1.96 hr (T VLM 513, lHallinan et all (12001 ). 3.07 hr (2MASS 
J0036+ 18. lHallinan et al.l (120081) ). and 2.83 hr (LSR 1835+32, 
lHallinan et al.l(l2008l) ). 

The L dwarf b inary, 2MASS 0746+20, reported by 
lAntonova et al.l (120081) from a mini-survey of UCDs at 4 . 9 GHz 
has a high mean flux level of 286+24 ^Jv. iBerger et al.l d2009l ) 
presented an 8.5 hr simultaneous radio, X-ray, UV, and optical 
observation of this binary. The strong radio emission consists 
mainly of short-duration periodic pulses at 4.86 GHz with P = 
124.32 ± 0.1 1 min. The radio pulses are 1/4 phase different from 
the Ho- emission. 

The narrow bunching of multiple pulses of both left- and 
right- 100% polarized radio emission detected from TVLM 513, 
which originate in regions of opposite magnetic polarity, reveal 
the likely pres ence of a dipolar c omponent to the large-scale 
magnetic field (Hal linan et al.ll2008l) . Zeeman Doppler Imaging 
(ZDI) observations have also shown that such a large-scale dipo- 
lar magnetic structure could exist on an M4 dwarf star, V374 
Peg, which is also a full y conyective rapid rotator similar to 
TVLM 513 (iDonati et al.ll2006t lMorin"etal]|2008h . The topol- 
ogy of magnetic fields on UCDs needs to be constrained by more 
observations. 

The radio emission composed of quiescent and pulsing com- 
ponents from UCDs can be associated with not only the ge- 
ometry of the emitting region and rotation of UCDs, but also 
the behavior of plasma in the magnetic field, i.e. the radia- 
tion mechanism. The very high brightness tempe rature and high 
(up to 100%) circular polarization of the pulses (Hallina n et al.l 
|2007|) point towards an efficient, coherent radiat ion mechanism, 
the el ectron cyclotron maser instability (ECMI, Melrose & Dulk 
(fl98l ). This wave magnification process of the free-space radi- 
ation modes could be induced by some kind of anisotropic ve- 
locity distribution s of electrons, such as a loss-cone distribution 
(lLau&Chull9 83). ring shell distribution, or horseshoe distribu- 
tion dPritchettlll984|) . 

ECMI was successf ully applied to the auroral kilometri c ra- 
diation (AKR) on Earth dWu & Led 1 9791: lErgun et al.l2000l) . de- 
cametric ra diation (DAM) on J upiter, Saturnian kilometric radia- 
tion (S KR) dZarkall9 98. 2004) and solar millisecond microwave 
spikes (Aschwanden 1990b). Various authors have suggested 
its presence in exoplanets dZarkall2007l: iGriefi meier et al.1 12007: 
Jardine & Cameron 2008]). The magnitude of any contribution 
from incoherent gyrosynchrotron or synchrotron radiation to the 
quiescent components in UCDs is uncertain. 



The generation of ECMI emission is dependent on the en- 
vironment of the emitting region, such as the magnetic field 
(strength, structure), the electron density distribution (number, 
energy), the line-of-sight source scale, and the angle betwee n 
the line of sight and the magnetic field (Mek-ose&Dulk][l982), 
and also the details (velocity spa ce gradients) of the loss cone 
distribution (Aschwanden 1990a). This results in the possibil- 
ity of transient radio emission. The findings of lAntonova et alj 
(2007) indicate that UCDs may also have sporadic long-term 
variability in their levels of quiescent radio activity. This phe- 
nomenon could be related to the change of the environment 
of the radio-emitting region, leading to self-quenching of the 
electron-cyclotron maser. Radio observations can help us to de- 
termine the magnetic configuration and the nature of the plasma 
in the emitting region. 

For the purpose of understanding the magnetic topology, we 
built an active region model based on the rotation of TVLM 513 
and the ECMI mechanism to simulate the observed light curve. 
We summarize the previous radio observation on TVLM 513 in 
§|2] In $3] we present the model, and results are given in $4] We 
discuss these results in $5] and draw conclusions in $61 



2. Previous observations on TVLM 513 

TVLM 513 is a young radio active M8.5V dwarf with a bolo- 
metric magnitu de of log(Lboi /L ) =» -3. 65, effective temper ature 
r eff « 2200K (lTinnev et al.l 1 993U 1 995tlLeggett et al.|2 001). and 
situated at a distance of d — 10.6 pc dDahn et al.l2 002). From the 
theory of the formation and evolution of UCDs and the absence 
of lithium on TVLM 513, it is reasonable to infer values for the 
mass and radius of this star of ~0.07 M and ~0.\R Q respectively 
dReid et alJl2002l:lChabrier & Baraffdl2000l) . 

With the VLA a highly right-circularly polarized (~65%) ra- 
dio event from TVLM 513 was detected with a flux density of 
~1 100 / / Jy, as well as persistent va riable emission at 8.46 GHz 
(Berger (2002)). Os ten et all d2006l) conducted a multifrequency 
VLA observation of TVLM 513 at 8.4, 4.8 and 1.4 GHz, using 
a strategy that involved time-sharing a single 10 hr observation 
between the various frequency bands. TVLM 513 was detected 
at each frequency band with only marginal confirmation of vari- 
ability and no detection of flares or strong circular polarization. 

Again, using the VLA, lHallinan et ail d2006h found per- 
sistent and periodic radio emission from TVLM 513 at 8.44 
GHz and 4.88 GHz simultaneously, with a period of ~2 hr. 
Subsequently, extremely regular periodic bursts (p = 1 .96 hr, up 
to ~4 mJy) of high brightness and high ly circularly po larized ra- 
dio emission were reported by Hallinan et al. (2007). Multiple 
bursts of both left and right 100% circularly polarized emission 
were detected. Interestingly, the radio emission can switch states 
from nearly 100% left polarization to 100% right polarization in 
each phase. 

Another r adio burst w i th a fl ux density up to ~4 mJy was 
presented by iBerger et al.l ([2008 ) from a period of simultane- 
ous radio, X-ray, ultraviolet, and optical spectroscopic observa- 
tions. Steady quiescent radio emission superposed with multi- 
ple, short-duration, highly polarized bursts was observed, but 
these authors reported a non-periodicity in the pulses/flaring ac- 
tivity. In a re -analysis of this dat a, plus data taken ^40 days later 
(June 2007). lDovle et all d2010l) reported the 1.96 hr periodicity 
in both datasets, deriving a more accurate period. 
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Fig. 1. A sketch of the active region model on ultracool dwarfs. 



3. Model 

In this section, we present a model (Figure [TJ to simulate the 
observational time profile of the radio flux density from TVLM 
513 by assuming ECM emission is the dominant radiation mech- 
anism. 

When plasma electrons are energized in magnetic flux tubes 
with converging legs and foot-points in a high density atmo- 
sphere (perhaps as a consequence of magnetic reconnection), 
some of these fast electrons collide with the high density at- 
mosphere and thermalise at the foot-points. The remaining fast 
electrons are reflected in the converging field by a magnetic mir- 
ror effect. This process results in the formation of the anisotropic 
distribution of the plasma in velocity space, i.e. a loss-cone dis- 
tribution. 

The plasma including a loss-cone distribution is unstable; in- 
stability arises very quickly from such a distribution. A large 
amount of free energy can be released via the insta bility and con- 
verted to electromagnetic waves - see iDulkl (119851) for a review 
of the process. An external magnetic-field-a ligned electric field 
dCattell et al.lll998b IZarkalll99ot lErgun et al.ll2000l) induc ed by a 
time-varying external current source rtOmura et al . 2003) would 
further modify the plasma velocity distribution to shell or horse- 
shoe form, leading to an enhanced ECMI emission. In this paper, 
we assume for simplicity that the ECMI emission from UCDs is 
driven by the loss-cone distribution. The existence of the electric 
field and its effect will be addressed in future work. 

In a loss-cone region where the ECMI operates, the maser 
radiation is c oncentrated on t h e surf ace of a hollow cone as 
discussed by Melr ose & DuUd d 19821) (see Fig. Q] top panel). 
The half-angle ar of the hollow cone depends on the ratio of 
the velocity of the plasma electrons to the speed of light, i.e. 
cosffo = v p /c. For example, if v p /c = 0.5, we have ao = 60°. 
The surface of the cone should be very thin with Aa « v p /c. 
The maser emission in the loss-cone region can be observed if 
the line of sight is located within a thin conical sheet with thick- 
ness Aa. For a low magnetic loop with a small angle between the 
magnetic field and the surface of the UCD, the maser emission 
can be seen when the emission is near the top of the loop. The 
maser emission from near the foot-point (where the magnetic 
field is almost perpendicular to the surface of UCD) can be seen 
when the loop is near the limb. However, for a large scale, the 
maser emission could have an angular distribution. In our cal- 



culation, we assume we can observe the maser emission in each 
flux tube in the active region. 

By analogy to the solar coronal radio emission powered by 
two populations of plasma from the Sun with different veloc- 
ity, we propose there be similar active regions on UCDs. We 
note that, however, an alternative mechanism, where the hot 
plasma beam could result from the interaction of a close-in com- 
panion of the UCD, i.e. magnetized or non-magnetized satel- 
lites resulting in auroral emission, similar to Io- Jupiter syst em 
(lOueinnec & Zarkall998trSaur et al.ll2004tlZarka et al1l2005l) . is 
possible. ECM emiss ion has been detected in c ompact objects, 
such as white dwarfs dWilles & Wull2004l2005l) or neutron stars 
(IWolszczan & F rail 1992]). We are not able to rule out this model 
and the competition between the two models should be investi- 
gated in further work. 

3.1. Flux calculation 

The fine structure and time interval of the observations described 
in $2] indicate that the radio-emitting region on TVLM 513 
would consist of complex substructures. The detected flux den- 
sity § may be the sum of several small sources 



Ml) 



(1) 



where Si is the flux density for a sm all source which can be 
determined by the relation dDulkll 1985b 



Hi! 



T h dQ. 



(2) 



where 7b is the brightness temperature of a source, k% is the 
Boltzmann constant, / is the observed frequency, c is the speed 
of light, d£2 is the differential solid angle. If we assume that the 
radiation is isotropic, the differential solid angle should depend 
on the radius of the flux tube r tu b e of the small source and its 
distance from the observer. The flux density S , can be expressed 
as 



Si = kv 



7b 



And 2 ' 



(3) 



For simplicity, we assume that the active region has a sym- 
metric shape, and the small sources are randomly distributed 
within the region. The number of small sources h(t) may be de- 
termined by rotation of the UCD. We set the time t — to be 
the moment when the first active region emerges in the field of 
view of the observer; this is also the onset time of a radio pulse. 
With the rotation of the UCD, the area of the active region seen 
by the observer increases until it reaches a maximum, and then it 
decreases and disappears from the field of view of the observer. 
The maximum number of small sources is 



At-2nR' cos 9/TucdY 



r tube 



(4) 



where Tucd is the rotation period, R' is the height of the active 
region, 9 is the latitude of the active region and At is the time 
interval from the beginning of a pulse to the maximum flux. The 
time duration of a pulse would be 2Af for a symmetric shape. We 
assume that the radio emission is from a thin shell of the active 
region near the surface of the UCD, so that we can take y « 2 
and R' » /?ucd where R\jcd is the radius of the UCD. 
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Fig. 2. The change of growth rate and brightness temperature with different plasma parameters and loss cone parameters. f p /f c = 
ratio of plasma frequency and gyrocyclotron frequency, T w = background wave energy, «h = density of hot plasma, 7h = temperature 
of hot plasma, a c = loss cone a ngle, N _ loss cone p itch angle distribution slope, 7b = brightness temperature. The range of each 
parameter is in the definition of lAschwandenl (Il990a). We plot each panel assuming 'standard' values of the other parameters (see 
TableQ]i. 



3.2. Brightness temperature and relative parameters 

The brightness temperature of small sources depends strongly 
on the growth rate of the ECM emission and incoherent radia- 
tion of the background plasma. If we assume the maser emis- 
sion is operated by the loss-cone distribution of a population of 
hot plasma expressed by a Maxwellian multiplied by the func- 
tion sin N (a/a c ■ n/2) and an isotropic Maxwellian distribution 
of a cold background plasma, the brightness temperature can be 



determined by the parameters of these two t ypes of plasma. We 
adopt the quasi-linear theory developed by Aschwanden ( 1990a) 
to determine the growth rate and the efficiency of energy conver- 
sion. Here, we summarize the basic assumptions and the effects 
of several free parameters. 

The quasi-linear code o f lAschwandenl (Tl 990b) describes the 
evolution of the ECM instability and the wave-particle inter- 
actions by solving the kinetic wave-particle equations in a lo- 
cally homogeneous plasma. The wave equation includes in- 
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Table 1. Parameters to determine the flux density. Plasma pa- 
rameters: «h = density of hot plasma, 7), = temperature of hot 
plasma, u = f p /f c = ratio of plasma frequency and gyrocyclotron 
frequency, n c = density of cold (background) plasma, T c = tem- 
perature of cold (background) plasma, T w = background wave 
energy; Loss cone parameters: a c = loss cone angle, N = loss 
cone pitch angle distribution slope, B = magnetic field strength; 
UCD parameters: R\jcd - radius of UCD, 7\jcd = rotation period 
of UCD, r tu be = radius of flux tube, 9 - latitude of radio-emitting 
region, A = size of emitting region. Note that the number density 
of cold plasma is associated with f p /f c . 
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duced gyroresonance emission/absorption (for the X-, O-, Z- 
magneto-ionic modes and whistlers), but neglects spontaneous 
emission, free-free absorption, collisional deflection, and spatial 
wave propagation. The coupled diffusion equation contains the 
quasi-linear diffusion coefficients due to maser growth/damping, 
but neglects slower processes like particle loss and source terms. 
This assumption corresponds to the strong diffusion case. The 
quasi-linear diffusion process of the ECM instability in the 
solar corona successfully accounted for the time profile and 
observational charac teristics of decimetric millisecond spikes 
(I Aschwandenl 1 990bh . 

In order to constrain the growth rate F, energy conversion 
factor e c and saturation of the maser, we need to know the initial 
nature of the two populations of plasma and the shape of the 
loss-cone, i.e. for the hot plasma: particle density n^, particle 
temperature 7V, for the cold plasma: particle density n c , particle 
temperature T c ; for the loss-cone: loss-cone angle a c , pitch angle 
distribution slope N. Rather than choowing values of the cold 
plasma density n c we fix its value via adoption of a value for 
u = fp/fc, the ratio of plasma frequency (/ p 9 x 10" 3 (n c )^ 2 
MHz) and gyrofrequency (/ c =s 2.86 x B MHz) where B is the 
magnetic field strength in G. 

In the quasi-linear diffusion process, the existence of back- 
ground electrom agnetic wave energy T w is taken into account. 
As discussed by lAschwanden ( 1990a) in the case of solar mil- 
lisecond spikes, the initial brightness temperature could be as 
low as the level of thermal bremsstrahlung in the range of 10 6 ~ 
10 8 K. However, because of gyroresonance or gyrosynchrotron 
radiation, an enhanced photon level could exist in a flaring loop 
which would affect the ECM process significantly, yielding a 
wave turbulence at the level of ~ 10 15 K. We apply T w in the 
range of 10 6 - 10 16 K in our model. 

The free energy of the plasma with a loss-cone distribution 
can be converted into an equivalent electromagnetic energy. The 
energy conversion factor e c can be defined as the ratio of the 
change of kinetic energy between the initial state and the final 
state of the plasma and the initial kinetic energy. In the present 
paper, the amount of converted energy is about 0.5%, i.e. the 
sa me as inlAschwandenldl990al) . 

I Aschwandenl d 1 990a) investigated how the parameters influ- 
ence the growth rate and brightness temperature generated by 



the ECM instability and gives the general formula for the key 
parameters as follows: 



T = 6.9x 10 



(l.25 x 10 6 cm- 3 )( 11 + 11 (30°) + (30°) 

l.5\-l 



x < 



°^(w) +o - 3o (i5y~ ) i 1 ' 15 ®- ' 15 ) 

^2.0 - 1 .2 (^j) + 0.2 ((fr) 2 | - 0.1 <m< 0.24 

0.19(1.0092- 0.0092 (^) 2 ), 0.24 < u < 1.0 

0.021(-1.416 + 0.440 ( ; f T )- 0.024 ((f;-) 2 ), 1.0 < u < 1.4 

(5) 



T h = 9.0 x 10 17 1 % 1''"' ( 1.4 - 0.4 f ^X 1 

1 1.25 x 10 6 cm- 3 / \ \30°/ 



,0.45 



. 1.2 



1.0 



logr, 



1.4-0.41 

0.05 



14 



(6) 



x < 



1. 



1400(£)- 3 . 



0.1 < u < 0.24 
0.24 < m < 1.0 
1.0 < u < 1.4 



We plot in Figure [2] the effect of various parameters on the 
growth rate and brightness temperature. 

4. Results and comparison with observations 

Many parameters, not only those describing the nature of the 
plasma but also those associated with the properties of the UCD, 
can affect the observable flux density significantly. In this sec- 
tion, we will discuss these parameters and compare our simu- 
lations with observations. In order to study the influence of the 
various physical quantities we adopt 'standard' values for each 
parameter, listed in TableQ] Unless otherwise noted, the simula- 
tions are always for magneto-ionic X-mode <x = -1 , harmonic 
number s — 1 . 

4.1. The effect of rotation 

In order to see the effect of rotation of the UCD, we first start 
a set of simulations by fixing the initial plasma and loss-cone 
parameters as listed in T able Q] In addition, as suggested by 
IChabrier & Baraffd (l2000l) , we adopt the radius (R VC d =0.1 R e ) 
as a constant in our simulations. We initially assume that the 
radio-emitting region is close to the equator (latitude = 30°). 

Given the above parameters, the size (A) of the emitting re- 
gion and the tube size of a small radio area, we show in Fig.[3]the 
influence of rotation. From the figure, we can see that the radio 
light curve is broadened with increasing rotation period (7\jcd) 
while the intensity does not change at all. This is because the 
intensity of the radio emission is only related to the behavior of 
plasma and the total size of the emitting region. 

Figure |4] shows the effect of the total size (A) of the radio- 
emitting region (left panel) and the radius (r tu b e ) of the flux tube 
(right panel). It is easy to understand that r tu b e cannot change the 
radio flux as much as 7ucd- It could, however, be a controlling 
factor for the observed oscillation of radio flux because each flux 
tube could have a different environment so that we may have to 
give a distribution for the initial plasma parameters (see 
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Fig. 3. The influence of rotation period (7\jcd) on the radio light 
curve. Vertical scale denotes the flux density while the horizontal 
scale indicates the time. Dashed line, solid line and dotted line 
are for 7\jcd 2.96, 1.96, 0.96 hr, respectively. 



A is an important parameter for the nature of the radio- 
emitting region on the UCD, since we can easily see from Fig. 
|4]that increasing its value can lead to a rise in both intensity and 
pulse duration. A can be constrained by the observations once 
we know the rotation period of the UCD. We are able to evaluate 
A from 

A ~ (At ■ IttRvcd cos 0/r UCD ) r (7) 

where the parameters are the same as in Eq. [4] For the standard 
parameters and y — 2, we find Af a 15 s. This implies that the 
observed Af can be used to estimate A. 

4.2. The effect of hot plasma in the loss-cone 

Because of the conservation of energy, the electromagnetic wave 
energy escaping from the radio- emitting region has come from 
the kinetic energy of the hot plasma. So the brightness tempera- 
ture of one flux tube should be proportional to the number den- 
sity and temperature of hot particles. 

We can see in panels (c) and (d) in Fig.[5]the variation of the 
flux density [S ) with the number density («h) and temperature 
(7h) of the hot plasma. Simulations in this section assume stan- 
dard values for all parameters except «h and T^. Increasing 7h 
from 10 7 to 10 8 K, S will increase by ~16 times, while S rises 
by a factor of 12 when «h is one order of magnitude higher. This 
is consistent with the change in the brightness temperature (see 
Fig.©. 

It is easy to understand that the population of hot electrons 
is shifted to higher velocities where the number of undamped 
resonance ellipses increases and thus the growth rate increases 
with increase in temperature. On the other hand the cold back- 
ground plasma will negate the loss cone and the growth rate 
will decrease sharply due to the lack of undamped resonance 
ellipses, if th e hot electron t emper ature approaches that of the 
cold plasma. iMelrose et alj ( 1 19841) gives a criterion for effec- 
tive cyclotron damping by background cold electrons, which is 



Th/T c < 10 ~ 20. This value was confirmed in the work of 
lAschwandenl(ll990al) . 

4.3. The effect of cold plasma and the initial wave energy 

The number density n c of the cold plasma can be obtained from 
the ratio of the plasma frequency to the cyclotron frequency 
(fp/fc)- Given a magnetic field strength (fi=1750 G), / c =4.9 
GHz. Since the standard value of/ p // c =0.1, / p =490 MHz. This 
implies n c = 3 x 10 9 cirT 3 , the standard value adopted in our sim- 
ulations. Note from Fig. [2] that changing f p /f c from say 0.1 to 
0.2 makes no difference to the brightness temperature. 

In panel (a) of Fig. [5] we show the effect of changing f p /f c 
for a constant B. We see that the flux density goes up a little 
from / p // c =0.23 to 0.5 (« c =1.6xl0 10 to 7.4xl0 10 cirT 3 ), and 
then drops when / p // c is increased further to 1.2 (« c =4.3xlO n 
cm -3 ). This result is consistent with the variation of the bright- 
ness temperature in Fig.|2]because the flux density in our simu- 
lation depends strongly on the brightness temperature. 

We note, however, that, although the growth rate declines 
persistently in Fig. [2] the brightness temperature increases 
steeply by almost one order of magnitude when the magneto- 
ionic mode changes from X-mode to (9-mode at harmonic num- 
ber s=l. Thi s is because of the Do ppler resonance condition as 
discussed by Aschwanden ( 1990a). The resonance ellipses cov- 
ering the unstable portion of the loss-cone distribution form a 
smaller region of positive growth for the (9-mode, and a higher 
wave level results from the same amount of energy conversion. 

We need to mention that only the dominant magneto-ionic 
modes, i.e. fundamental (s=l) X-mode when / p // c < 0.24, fun- 
damental (s-l) (9-mode when 0.24 < / p // c < 1.0, and second 
harmonic (s=2) X-mode when 1.0 < / p // c < 1.4, have been 
taken into account in our simulations. When f v /f c > 1.4, some 
electrostatic instabilities become important rather than the ECM 
instability, and a condition for ECM emission escaping from 
plasma is / p // c «C 1. The fundamental (s-l) Z- mode may be 
domi n ant when f D / fc ~ -3 and fp/f c ~ 1.1 (IMelrose et al.l 
11984 lAschwandenT ll990al) . However, we omit this mode be- 
cause its appearance depends strongly on the ratio / p // c and an- 
other mechanism would be needed to convert it to electromag- 
netic radiation. 

When T w (defined as the initial wave energy) results mainly 
from thermal bremsstrahlung it takes a low value of 10 6 ~ 10 s 
K. J w may become larger, in the range of 10 10 ~ 10 16 K, due to 
incoherent gyroresonance or gyrosynchrotron radiation. Lower 
T w can lead to a higher energy conversion efficiency because 
of the effect on the resonance ellipses. The flux density drops 
slightly with increasing T w as seen in panel (b) of Fig. [5] 

4.4. The effect of loss-cone parameters 

Panel (e) in Fig. [5] shows the effect of the loss-cone angle on 
flux density. The rapid increase of flux density when the loss- 
cone angle (a c ) changes from 10° to 20° is the consequence 
of the significant increase in the size of the region of velocity 
space involved. Quasi-linear diffusion influences a large num- 
ber of particles and results in growth and amplification of some 
wave modes growing. On the other hand, when a c > 20°, the 
positive gradient in the perpendicular direction decreases so that 
the electromagnetic wave energy density can not be amplified 
effectively. 

The effect of pitch angle distribution slope is similar to that 
of the loss-cone angle, see panel (f) in Fig. [5] 
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Fig. 4. Flux density against time with different size of the emission region (A, left panel) and radius (r tu b e , right panel) of one flux 
tube. 



4.5. Determining the environment of the radio-emitting region 

In this section, we determine the possible range of the different 
parameters. From Table [T] there are at least 13 parameters that 
can influence the radio light curve of TVLM 513; we are able to 
reduce this number, however. 

Considering the contribution of the degenerate electron 
gas and the ionic Coulomb pressure, the radius is almost 
a constant around 0.1 R G , in the range of O.O8/? to 0.11/? Q 
dChabrier & Baraffell2000l) . 

The temperature of cold plasma is more difficult to deter- 
mine. In this work we assume it is around 10 6 K, similar to the 
typical temperature of the solar corona. On the other hand, if 
the coronal temperature was lower, this would favor wave prop- 
agation rather than damping. We set the initial wave energy T w 
at a high level of 10 14 K. On the other hand, from Fig. \5\ the 
flux density varies only slightly (a factor of 2) wh en T w changes 
from 10 10 to 10 15 K. We should however note that lHallinan et al.1 
(2008) suggested that the ECM instability may be a viable source 
of quiescent unpolarized radio emission, indistinguishable in 
temporal and polarization characteristics from gyrosynchrotron 
radiation. 

For convenience we adopt a plausible value of 0.1 for the 
ratio u of plasma frequency to gyrofrequency (/ p // c ), a value 
that permits unstable growth of the ECMI. If we focus on the 
observation frequency of 4.9 GHz, the magnetic field strength 
{B) has to be around 1750 G and the density of the cold back- 
ground plasma is ~3xl0 9 cirT 3 . Any change of u alone will not 
affect the maximum flux density for the fundamental X-mode 
and fundamental (9-mode respectively. However, we stress that 
if we change the magnetic field strength (B), even keeping u at 
the same value, the maximum flux density will change because 
it depends on the cyclotron frequency. 

The size of the radio-emitting region A is associated with the 
time duration of the radio pulses. Thus in the case of TVLM 513, 
we have A * (54.75Af) r when the latitude of the active region 



is 30°. Changing the latitude to 70° gives A * (21.62Af) r . At is 
easily obtained from observations and is a more convenient pa- 
rameter than A. r tu b e is a variable to describe the radius of a flux 
tube. We expect there are distributions for plasma parameters 
(Th, «h) and loss-cone parameters (a c , N) in different flux tubes. 
The variation of different parameters can result in the oscilla- 
tion of the observed flux density. Smaller values for r tu b e lead to 
smoother simulated radio light curves (see Fig. |4j). 

So now the free parameters are reduced to 6, i.e. the density 
«h and temperature 7h of hot plasma, the angle a c and pitch angle 
distribution slope N of the loss-cone, the radius r tu b e of one flux 
tube and the size of the emitting region A. The functions of the 
parameters are: «h and Th can control the change of flux density 
dramatically, a c and N can control the change of flux density 
gently, A can control both flux density and time duration of the 
radio pulse which can be constrained by observations, and r tu b e 
describes the radius of the flux tube in the active region. 

4.6. Comparison with observations 

The free parameters can be constrained by making our simu- 
lations match up to observations. Fig. [6] shows the comparison 
of our simulation with observations (see Hallinan et al. 2010 
for further details). The observation is taken from 2008 May 
19 (UT) at 4725 MHz for the M8.5V dwarf TVLM 513. This 
means that the magnetic field strength is 1652 G because of the 
cyclotron frequency. We adopt f p /f c =0.1, implying the density 
of cold plasma in the emitting region is ~2.64xl0 9 cirT 3 . The 
size of the radio-emitting region is constrained by the time du- 
ration of the two pulses (A w (54.75Af) 2 ), implying 1620x1620 
km 2 for the first pulse and 2060x2060 km 2 for the second pulse. 
The radius of each flux tube is taken as 180 km arbitrarily be- 
cause we are not able to get any information on its value from 
the observed radio light curve. We take a c = 30° and N — 5 for 
the first pulse, a c — 35° and N — 2 for the second pulse. 7h and 
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Fig. 5. Simulated maximum flux density of a radio pulse as a function of the ratio of cyclotron frequency and plasma frequency 
(/p//o)i background wave energy (T w ), density of hot plasma («hX temperature of hot plasma (7h), loss cone angle (a c ), loss cone 
pitch angle distribution slope (AO are shown in panel (a) to (f) respectively. The parameters for the calculation in this plot are the 
same as in Fig. [2] 



«h are same for the both pulses, equal to 10 7 K and 1.25 x 10 5 
cm 4 respectively. 

The simulated light curve reproduces the observation. 
However, this does not mean that the above values for the differ- 
ent parameters are unique. In order to compare our simulation 
with the observed two pulses, we take a flat distribution for a c 
and N in the range of 10° - 50° and 1-6 for the both pulses, 
n h : (1.25 - 5) x 10 s cirT 3 , T b : (1 - 4.5) x 10 7 K for the first 
pulse and « h : (1.25 - 4.5) x 10 5 cirT 3 , T b : (1 - 3) X 10 7 K 
for the second pulse (see inset panel in Fig . [6j> . We note that the 
size of the emitting region depends on the latitude 9 in our sim- 
ulation (see 34.5b . Changing the latitude from 30° (the value for 
the simulation) to 70° gives a value for the size of 640x640 km 2 
for the first pulse and 814x814 km 2 for the second pulse. Other 
parameters also need to change in order to fit the observed light 
curve. 

The simulation in Fig. [6] is remarkably similar to the obser- 
vations. Our results indicate that the two radio emission regions 
producing the two pulses are very close, which would be consis- 
tent with the nature of the cone radiation of ECMI. An important 
feature is that the radio pulses would repeat with the rotation pe- 
riod of the UCD. In the case of TVLM 513, the period is -1.96 
hrs. Moreover, we note that the decay time of the pulses in some 
observations is longer than our simulation. This is possibly due 
to the deformation of the radio-emitting region as fast rotation 



of the dwarf can cause the shape of the emitting region to vary 
from an almost symmetric circle to an asymmetric ellipse with a 
tail. 



5. Discussion 

Our results show that rotation coupled with the ECMI mech- 
anism can account for the flux density and polarization of the 
radio pulses from TVLM 513 successfully. We can not exclude 
the possibility that the depolarization of ECMI could be due to 
radiation transfer of the emission in a neutral atmosphere with 
lower fractional ion ization or that inhomogeneous dust clouds 
( Littl efair et al.l2008h could have an eff ect on the quiescent emis- 
sion and the unpolarized components. lHallinan et al.l |2006) sug- 
gested that the depolarization or mode conversion of the X-mode 
emission occurs in a density cavity, as mode conversion of ter- 
restrial kilometric radiation (TKR) from X-mode to R-mode in 
the emitting density cavity (Ergun et al. 2000) may account for 
escape of maser emission without re-absorption at higher har- 
monics of the emission frequency. 

For the optically thick source, i.e. optical depths r » 1, the 
brightness temperature T\, can be constrained by Eq.[3]by using 
the observed flux density. In order to see the relation between T\, 
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Fig. 6. Comparison of our simulation with observation taken with the Arecibo telescope by Hallinan et al. (2010). Black solid line 
shows the observation of 2008 May 19 (UT) at 4775 MHz. Red line shows the simulated curve. The inset panel shows the simulated 
light curve using a uniform distribution of the free parameters. 



and flux density, we rewrite the equation in the form of 



T h = 5.32 x 10 10 K 



f 



4 mJy / \4.9 GHz/ \10pc 



Rs_ 
Rjup 



(8) 

where S / is the flux density in mJy at the frequency / (GHz), d 
the distance of the radio source from us in pc, /? s the size of the 
emission region in Jup i ter radii (1 j?j UD ~ 0.1 fi R »7x 10 9 cm) 
dLinskv & Garvl [l98l iDovle et all U98 8I: iDorman etaf] [19891 
Burrows et al. 1989t iLeto et al.ll2.000h . Unfortunately, when we 
calculate 7b, we have to assume the size of the radio-emitting re- 
gion. For example, observations show that the flu x density of the 
pulses of TVLM 513 is about4 mJy at -4.9 GHz. lBergerf [2002) 
obtained the brightness temperature in the range of 10 8 - 10 9 
K by assuming the size of a corona to be 2 - 4 Rj up , while 
lHallinan et all (120061) deduced a value of 2.9xl0 10 K if the size 
of the region is 1 Rj up . In our model, It, > 5 x 10 10 K is about the 
temperature of the quiescent radio emission, assumed to come 
from a large emission region (< 1 Rj up ). For the radio pulses 
from magnetic loops, the theoretical temperature of the coher- 



ent ECMI emission can be up to 10 15 K, implying the emission 
region for the pulses is much more compact, e.g. 0.007 Rj up . 

The configuration and topology of magnetic field on UCDs 
remains unclear. In our model, we only need a simple dipole 
poloidal-like field to calculate the flux density and explain the 
high polarisation of the radio pulses. In the case of TVLM 513, 
multiple bursts of both left and right 100% circularly polarized 
emission in regions of opposite magnetic polari ty indicate the 
existe nce of a dipolar large-scale magnetic field (Hallinan et al. 
[2007h or a few small active regions with scale << 1 Rj up . 

Another important parameter, when combined with the mag- 
netic field strength, is the pitch angle /3. Electrons with small f3 
less than a critical value /3 C , i.e. /3 < f3 c , precipitating into the 
dense atmosphere are lost, while the electrons with ft > /3 C will 
be reflected back to the flux tube to form an anisotropic velocity 
distribution. The value of f3 c depends on the convergence factor 
which is determined by the ratio of magnetic field strengths at 
the top and the foot-points of the flux tube (also c alled magne tic 
mirror ratio). For a symmetric flux tube, we have 

dDuMl985) 



/3 C = arcsin(B top /B foot ) I/2 



(9) 
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Fig. 7. Possible magnetic field structures for UCDs. The solid 
line is for Eq.fTUl which was suggested for a dM star, the dashed 
line for a unipolar-like field, the dotted line for the possible field 
in the solar coronal region. 

where B top and Bf oot are the magnetic field strengths at the top 
and in the foot-point of the magnetic flux tube, respectively. 
Typical values of B top /Bf oot in the Sun are in the range of 0.1 to 
0.5. In the simulation, the pitch angle is approximately equal to 
30°, which means we should have 0.5 > sinyS c = {B top /B{ oot ) 1 ^ 2 . 
This could give a lower limit on Bf oot . 

Assuming the magnetic field is radial, the field strength of a 
region can be described by 

B = T(Hf • 

where Bq is the magnetic field at the photosphere, R the height 
of the radio emitting region measured from the centre of the star 
and R t the photospheric radius. This was chosen as a compro- 
mise between a unipolar-like field B = Bo(R/R*)~ 2 and the best- 
fit form for solar coronal magne tic fields above active regions 
B = 0.5 [(R/R,) -I]' 15 given bv lPufk & McLean! (1 19781) for the 
range 1.01 < R/R* < 10. This choice makes the field fall off 
more quickly than either a unipolar magnetic fie ld configuration 
or tha t of the solar active region magnetic field (Gary & Linskvl 
11981b . Figure|7]shows the dependence of magnetic field strength 
on height of the radio emission region above the photosphere. 

Combining Eqs. |9] and [TUJ the radio emission region and 
therefore the magnetosphere can be constrained in the range of 
0.56 R* to 6.75 R*. A possible position for the radio-emitting re- 
gion in the radio active UCDs at ~4.9 GHz is 1 .5 R, from their 
center, 1.06xl0 5 km in the case of TVLM 513. This implies that 
the magnetic field strength at the photosphere or chromosphere 
would be as large as 7,000 G. More observations in the optical 
(Zeeman Doppler effect) and infrared bands are needed to con- 
strain the magnetic topology and the behavior of the plasma in 
the magnetic flux tube and to determine whether such large field 
strengths exist on these objects. 

In addition, we expect the existence of an enhanced ambi- 
ent wave energy background by gyroresonance turbulence or gy- 
rosynchoroton radiation and some intense events at other wave- 
lengths, e.g. optical or X-ray emission, which would occur from 
the process where hot plasma starting from collision-less region 



collides with the collisional chromosphere or photosphere. In the 
case of X-ray emission, thermal bremsstrahlung emission and 
inverse Compton scattering could be the responsible mechanism 
since there are hot plasmas with 7h up to 10 8 K and possible low 
energy photons. However, since the Thomson scattering optical 
depth of the corona of the brown dwarf is <SC 1, the contribution 
of inverse Compton scattering is unlikely to be important. 

For the thermal bremsstrahlung emission of ionized hy- 
drogen and helium dominated source, the detected flux 
density integrated over frequency is Sx-b = 4.67 x 
lO~ 2S T^ 2 n^n c Z 2m g B Rld^ 2 , where Z is an ion of charge in units 
e (here we take Z = 1), g B the Gaunt factor (we take g B = 1.2, 
which gives an accuracy of <20% since 1.1 < g B < 1.5), R s 
the radius of the source, d d istance of the source from us and 
all quantities are in cgs units (Rvb icki & Lightmanll979l) . In the 
case of TVLM 513, assuming that the X-ray emission comes 
from the same region as the radio emission, we can take the 
parameters as T h = 10 7 K, n h = 10 8 cm 4 , n c = 3 X 10 9 
cm" 3 , R, = 2 x 10 8 cm, d - 10 pcx 3.1 x 10 19 cm, hence 
we get the Sx-b = 4.42 x 10~ 21 erg-cm _2 -s _1 , giving an X- 
ray luminosity of Lx-b ~ 5.34 x 10 19 erg-s -1 . The X-ray flux 
density/luminosity could be underestimated significantly as the 
X-ray emission might be diffuse and from a larger region (per- 
haps -10 times) than that of the radio emission, as is t he case for 
the X -ray emission observed from Jupiter by S uzaku (lEzoe et all 
l2010l) . 

Interestingly, BergeretaTJ {2008) obtained a marginal de- 
tection in X-rays suggesting a flux density of 6.3 x 10~ 16 
erg-cirT 2 ^ 1 (luminosity Lx - 8.5 x 10 24 erg-s~') with mean en- 
ergy at 0.9 keV. This means the temperature of the hot plasma 
could be slightly lower than 10 7 K. Further multi-wavelength 
observations will help to refine our model and its parameters to 
understand the radio and X-ray emission from these kinds of 
cool objects. 



6. Conclusions 

An active region model is applied to the radio emission from a 
cool dwarf, in which the ECMI mechanism is responsible for 
the radio bursts from the magnetic tubes, while the rotation of 
the dwarf can modulate the total observed flux with respect to 
time. The time profile of the radio light curve is in the form of 
power law in our model. Using this model, we can determine 
the nature (e.g. size, temperature, density) of the radio-emitting 
region plus the magnetic topology can be constrained as well. 

In the case of TVLM 513, our model shows the loss-cone 
electrons have a density in the range of 1.25 x 10 5 - 5 x 10 5 
cm" 3 and temperature between 10 7 and 5 x 10 7 K. The brightness 
temperature is typically ~ 10 15 K for pulses, ~ 5 x 10 10 K for the 
background emission, implying the ECMI mechanism operates 
in compact region of ~0.007 Rj up if the active region is at 30°. 
For an active region closer to the pole, e.g. 70°, the size is ~60% 
smaller, implying a higher brightness temperature. 

The model predicts an enhanced ambient wave energy back- 
ground and a ^7000 G surface magnetic field strength. The the- 
oretical X-ray flux density in our model is much smaller than 
a marginal X-ray observation of TVLM 513, which implies a 
more complicated plasma behavior or magnetic structure on the 
dwarf. Additional multi-wavelength observations are needed to 
constrain the tentative conclusions and help us to improve the 
understanding of the magnetic field on ultracool dwarfs and to 
test the viability of this model in comparison with others, such 
as the auroral model. 
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